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Abstract 

It is found that the seesaw mechanism not only explain the smallness of neutrino masses but also 
account for the large mixing angles simultaneously, even if the unification of the neutrino Dirac 
mass matrix with that of up-type quark sector is realized. We show that provided the Majorana 
masses have hierarchical structure as is seen in the up-type quark sector and all mass matrices 
are real, we can reduce the information about the absolute values of neutrino masses through the 
data set of neutrino experiments. Especially for ^13 = 0, we found that the neutrino masses are 
decided as mi : 7712 : mj, ~ 1 : 3 : 17 or 1 : 50 : 250 {mi ~ m2 ■ ^3 ~ 3 : 1 or 12 : 1) in the 
case of normal mass spectrum (inverted mass spectrum), and the greatest Majorana mass turns 
out to be m§ = 1 x 10^^ GeV which just corresponds to the GUT scale. Including the decoupling 
effects caused by three singlet neutrinos, we also perform a renormalization group analysis to fix 
the neutrino Yukawa coupling matrix at low energy. 

PACS numbers: 14.60.Pq, 14.60.St, 12.15.Ff 
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I. INTRODUCTION 



Neutrino sector has many curious properties which are not shared the quark and 
charged lepton sectors. For example, neutrino masses are very small [1| compared with 
those of quarks and charge leptons. The large mixing angles seen in the experiments of 
atmospheric neutrino oscillation and long baseline reactor neutrino oscillation (related to 
solar neutrino deficit) y, H are also the new feature, not seen in the quark sector. 

nnn 

It is well known that the seesaw mechanism p, [71, |8[ can explain the smallness of neutrino 
masses naturally. In this mechanism, neutrino mass matrix which describes the low energy 
observables is given approximately by 



ID, 



fl.ll 



where Aio and Ain are Dirac and Majorana mass matrices of neutrino, respectively. If we 
require that the order of magnitude of A^d is the weak scale and that of A4r is the GUT 
scale, we can roughly obtain the desired order of magnitude of Ai^- 

In addition, it was pointed out in some articles |9| that this mechanism also may be 
responsible for the enhancement of the mixings in the leptonic sector, compared with those 
in quark sector. This enhancement mechanism can be seen in the case of simplified two- 
generation scheme as follows. If we assume that is symmetric and real matrix and Mr 
is real matrix (and necessarily symmetric), we can represent A4d and Ai^^ as 
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where r = mf/m^, R = mf/m^ and ce = cos 6, se = sin 6 etc. We have assumed that 
Ai D coincides with that of up-type quarks in the basis where down-type quarks are in their 
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mass eigenstate. Namely 6 is the Cabibbo angle. M.,y is then represented as 
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This means that in the leptonic sector the mixing angle is replaced as 6 
condition 

R-r^ + tan^(6' - 0) < r tan(6' - 0) 



(1.5) 

6 -\- a and if the 
(1.6) 



is satisfied, a becomes large considerably and the enhancement of mixing angle is realized 
due to the seesaw mechanism. In this context, the large mixing angle 6' + a is induced only 
for certain relation between three parameters, i.e. r, i?, tan(^ — 0). 

Inspired by this observation, in what follows, we apply this mechanism to the case of 
realistic three-generation scheme, assuming that is a symmetric matrix and all mass 
matrices are real. The former assumption is justified in the case of 50(10) GUT, because 
this model contains a subgroup SU{2)l x SU{2)r x U{1)b-l and has left-right symmetry. 
Since the exchange of chirality L ^ R corresponds io M.d ^ -^dj left-right symmetry 
implies that A^^i is a symmetric matrix, if we neglect the radiative corrections completely. 
The framework of 5*0(10) GUT, we work on, will be discussed in some detail in Sec lIVI The 
l.h.s. of Eq. ljl.lll . being a real symmetric matrix, has 6 degrees of freedom which are related 
to the low energy observables, i.e. 3 mixing angles and 3 mass eigenvalues. On the other 
hand, the r.h.s. of eq ljl.H) is parametrized by 12 parameters, i.e. 6 (3+3) mixing angles 
and 6 (3+3) mass eigenvalues of Dirac and Majorana mass matrices. These mean that the 
observables are not enough to determine the parameters. If, however, we regard A^d as 
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that of up-type quark sector, as implied by 5'O(10) GUT, and assume that the mass ratios 
of Majorana masses have a hierarchical structure in analogy to what we observe in quark 
sector, i.e. 

mf : : = : R:l, (1.7) 

the number of parameters is reduced to 5, so that the 5 parameters can be fixed by 5 
observables in neutrino oscillation. Namely, if we use the following 5 observables |lfl| . 

tan'Osoi = 0.39^°:°^ (1.8) 
tan^ Oatn. = l+ai (1-9) 
sin^^is < 0.041 (1.10) 
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X 10^^ eV^ (1.11) 
X 10"^ eV^ (1.12) 



the 5 parameters are fixed and we can finally estimate a remaining observable, i.e. the 
absolute value of neutrino mass. Identification oi M.^ with that of up-type quark sector is 
justified only at the GUT scale, not at low energies. We perform a renormalization group 
(RG) improved analysis in Sec llVl 

This paper is organized as follows. In SecjUl we parametrize the mass matrices in the 
realistic three-generation scheme, and discuss how the mixing parameters out of 5 parameters 
are describable in terms of R with the help of neutrino oscillation data. In Sec lIIIl we utilize 
the remaining datum of neutrino oscillation to determine R and finally derive the absolute 
values of neutrino masses. In Sec lIVl we improve these estimations taking into account the 
renormalization group (RG) effects to M-o including the effect of successive decoupling of 
three heavy singlet neutrinos. SecElis devoted to the conclusion. 

II. THE THREE-GENERATION SCHEME 

In this section, we discuss the realistic three-generation scheme. We adopt the following 
minimal set of assumptions: 

• the Dirac mass matrix, M-d/vs, symmetric and real. 

• the Majorana mass matrix, TW/j, is real. 



Ai D coincides with the mass matrix of up-type quark sector at the weak scale. 



• the eigenvalues of Majorana mass matrix have a hierarchical structure and are well 
approximated by mf : mf^ : m§ = B? : R : 1. 

In what follows, we adopt the basis in which down-type quarks and charged leptons are 



in their mass eigenstates, as suggested by 5*0(10) GUT, and use the fol 
Dirac mass matrix of neutrino at the weak scale (the 3rd assumption) 
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A. Parametrization 



According to the above assumptions, Dirac mass matrix is parametrized as 
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where we use an approximation that the Dirac neutrino masses are represented as 
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TJlT : 7712 • "^3 ~ 



(1/290)^ : 1/290 : 1 = : r : 1. 



As an explicit representation of v, we use a standard parametrization |12l.ll3. 
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where cu stands for cos^fg ^^c. 

On the other hand, Majorana mass matrix is represented as 



u 



1 



U 



1 



u 



R 



u 



(2.6) 



where u is an orthogonal matrix. 

Combining Eq. ()2.3|) and Eq. ()2.6|) . we can represent neutrino mass matrix as 
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where V = u^v denotes the deviation of u from CKM matrix, and contains 3 mixing angles. 
On the other hand, l.h.s of Eq. (|2.7|) can be written by observables as 



M, = U 



mi 
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(2.8) 



where mj's are neutrino masses and U is the MNS matrix except for the CP phase. 



U = U23U13UU 
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(2.9) 



cos^^i3 sin^^i3 

1 
- sin 6*13 cos ^^13 

These observables have already got certain values or constraints as seen in 

Eas. finni(fni|i(fTTn|)(im|i(frT^ . 

For a while we concentrate on the mixing angles and mass ratios, ignoring (m|^)^/mf in 
Eq. (|2.7|) . Then the following 4 observables obtained from the neutrino oscillation data 



tan^6' 



12 



0.39 



(2.10) 



23 



13 



ms — m 



irrin 



mil 





8.2 X 10^5 
2.2 X 10-3^ 



(2.11) 
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are enough to fix the (3 angles in the) mixing matrix V and the mass ratio R. 



B. Constraints to mixing matrix 

To decide the mixing matrix V and R, we define a matrix A by taking off a factor 
(mf )2/mf from A^,, in Eq. (jT7|) 
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All Ai2 Ai3 
Al2 A22 A23 
Ai3 A23 A33 



(2.15) 



As this matrix is proportional to Aiu, it should be diagonalized by MNS matrix U. Namely, 



A (X U 



mi 



m,2 



1 



(2.16) 



Then by using Eq.^JU\) to Eq.^J^, i.e. 623 = 45°, ^13 = 0° and 6^ = 6 = 32° (this 
6 should not be confused with that used in the two-generation scheme), for the mixing in 
Eq. ()2.9jl . we get from this relation 3 conditions to be satisfied by the elements of A: 
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We also have relations to express mass ratios, pi, P2, in terms of the elements of A: 
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Basically, the 3 conditions in Eqs. ()2.17j) . ()2.18|l . ()2.19|l can be used to express the 3 mixing 
angles in V, i.e. 0235 013,012 in = V23V13V12, in terms of R. Then the matrix A is 
completely written in terms of i?, and substituting Eqs. ()2.20|) . 1)2.211) in Eq. ()2.13|) . R can 
be fixed. Using Eq. ()l.ll|)()1.12|) . we can finally determine all absolute values of neutrino 
masses, mi, 1712 and m^. The matrix A, however, has a complicated form when written in 
terms of (pu, 023, 0i3 and R, so that it is not easy to represent V in terms of R. To avoid 
this difficulty, we will utilize the orthonormality of V as much as possible, and will use only 
a minimal approximation. 

C. Appropriate approximation 

It is worth noting that, because of the orthonormality of V, for R = 1 the matrix V 
becomes irrelevant, since 



where A = diag{l, R, R'^). To achieve a minimal approximation mentioned above, let us 
define a matrix X as 



and discuss about the relations between X and R, instead of discussing 0i2, 023,013 them- 
selves. As A is a symmetric matrix described by V2i, (i = 1 — 3), it contains 6 components, 
which we denote as Ai to Aq: Ai = An, A2 = A12, A3 = A13, A4 = A22, A5 = A23, Ag = 
A33. At first glance, one may think that the degrees of freedom in both sides of Eq. ()2.23)l 
don't coincide. We, however, can avoid this mismatch utilizing the orthonormality of V as 
is seen below. 

Using above definition and Eq. ()2.14|) . each element of A may be written as 



V^AV = V^{I -{I - X))V = 1- V^{I - A)V, 
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A = V^{I - A)V, 



(2.23) 




(2.24) 



k=l 



8 



where dij S cLFG defined as 
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Then the three conditions, Eqs. ()2.17j) . ()2.18|l . ()2.19|l . lead to coupled linear equations for 

. A22 = A33 - 4?)Xfc = 4? - (2.32) 

k=l 

. A,2 = -A,3 ^ E(4? + a['3^)Xk = aS? + aj? (2.33) 

k=l 

tan 20 IZ(an - 4? + ^2? + 2^2 cot 29a^S)Xk 
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= afi - ^22 + ^2? + 2 V2 cot 2eafj . (2.34) 
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As mentioned above, not all of Xi are independent by orthonormality. For instance, due 
to the orthonormality of V we find Tr(X) = Tr(/ — A) = l — — i?^, i.e. 

Xi + X4 + Xe = 1 - + 1 - i?^ (2.35) 

Similarly Tr(X^) yields 

Xl ^2X1 ^2X1^X1^ 2X1^X1 = {1- + {1- RY- (2.36) 

Another relation comes from the determinant of X: 

XiX^Xq — XiX^ — X4X3 — X5X2 ~l~ 2X2X3X5 = 0. (2.37) 

In this way, we have found, in total, 6 coupled equations for Xj, which are enough to solve 
for Xj in terms of R. However, two of them are non-linear equations and are very difficult 
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to solve in terms of R analytically js^. Alternatively, we will make one approximation in 
Eq. ()2.35j) and divide it into two equations. Using Eq. ()2.23|) . we can express Xi + X4 as 

+ X4 = (1 - R) {1 + (1 + R)V^3 + RVii) ■ 

Assuming that 1^13 -C 1 and -C 1, we use following approximated expressions. 

Xi + Xi ~ l-R (2.38) 
Xg ~ (2.39) 

When these equations are combined with the linear three conditions (j2 . 32|) . (j2 . 33|) . (j2 . 34p . we 
can represent X2,X3,X4,X5 and Xq in terms of linear combinations of Xi and R. Sub- 
stituting these in the remaining non-linear equation ()2.36|) . we obtain a quadratic equation 
about Xi which can be solved easily and finally obtain the expressions for Xi, ■ ■ ■ ,Xq as 
functions of R. We can finally fix the value of R by utilizing the remaining condition ()2.13p . 
i.e. 

pI{r)-pI{r) ^ 8.2 X 10-5 
|i -pi(i?)| 2.2 X 10-3' ^ ■ ' 

Note that there are two possible cases reflecting the uncertainty of the sign of mass-squared 
difference in atmospheric neutrino oscillation experiment, i.e. normal or inverted mass spec- 
trum, i.e. Am^^^ > or Am^^^ < 0, respectively. There are some proposals to fix the sign 
of atmospheric neutrino mass squared difference, i.e. discrimination between normal mass 
spectrum and inverted one by utilizing the difference of matter effect of the earth between 
electron neutrino and electron anti-neutrino at Neutrino Factorvjl?!. llS^. 

III. NUMERICAL RESULTS 
A. Normal mass spectrum 

In the case of normal mass spectrum, two values of R are allowed 

R = (3.214, -3.780) x 10"^ = (2.703, -3.179) x (3.1) 

Substituting these values in Eqs. (j2.20j) . (j2.2ip . we can immediately calculate the values of pi 
and p2, i.e. the ratios of neutrino masses 

= = (_0.003977, -0.05856), (3.2) 
10 



P2 = — = (0.1896,0.1981). (3.3) 

1713 

We can interpret the sign of mass eigenvalues seen in Eqs. ()3.2|l . ()3.3|l including the overall 
one in Eq. ()2.7|) as CP-parity in neutrino sector and can define all of mass eigenvalues as 
positive without a confiict with CP-invariance. For detailed discussion about CP-invariance 

a„d CP-pant, see .efs.QQ. 

Taking into account this fact, we can eventually determine the absolute values of neutrino 
masses, invoking to the best fit values of Eqs. (jl.ll|) . (11.121) : 

mi = (0.1900,2.802) x 10~^ eV (3.4) 
ma = (9.057,9.479) x 10"=^ eV (3.5) 
ms = (47.77,47.85) x 10"^ eV. (3.6) 

Furthermore, from Eqs. (j2.7|) . (j2.8|) and mf = 181 GeV, we can estimate the greatest Majo- 
rana mass as 

m3 it^ 

Substituting the obtained values of R to Eq. ()2.37p . we can immediately calculate the 
values of det{X) and we can find that it does not vanish exactly but has an extremely small 
value, O{10~^^), which can be regarded as the error of the numerical calculation in the last 
step to fix R. We also can confirm that the approximate relation, ()2.38j) . ()2.39|l . hold within 
this error. 

It is worth mentioning that the absolute values of neutrino masses at low energy obtained 
above are similar to the results obtained in Ref. 0|, in which the analysis are based on 
SUSY 5*0(10) X SU{2)p. They make use of quark-lepton symmetry and left-right symmetry 
at Mgut relying on 5*0(10), and a specific mass matrices with five texture zeros due to 
SU{2)p family symmetry. Though we also make use of quark-lepton symmetry and left- 
right symmetry, we do not utilize the family symmetry and assume that the Majonara 
masses have a hierarchical structure ("geometric mass hierarchy" j2^) instead. 
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B. Inverted mass spectrum 



We can just follow the procedure given above for the normal mass spectrum, and find 
two values of R 

R = (2.927, 6.027) x 10"^ = (0.2462, 0.5069) x (3.8) 
The values of pi and p2 take 

= = (11.62,2.739), (3.9) 

P2 = — = (11.84,2.785), (3.10) 

and we eventually get the absolute values of neutrino masses, 

mi = (4.619,4.943) x 10"^ eV (3.11) 

m2 = (4.707,5.026) x 10"^ eV (3.12) 

mg = (0.3976,1.805) x 10"^ eV, (3.13) 

and the corresponding the greatest Majorana mass 

= (14.70, 3.861) x 10^^ GeV. (3.14) 



C. In the case of sin^ 013 < 0.041 

For simplicity, we have used the assumption 613 = in Eq. (j2.12p . as ^13 is stringently 
constrained. It, however, will be desirable to evaluate the dependence of mass eigenvalues 
on ^^13 at least for its 3cr allowed range, i.e. sin^ ^13 < 0.041. There is no difficulty to achieve 
this, once we employ the same prescription as the one mentioned above. We found that 
generally the solution for mi 2,3 are multi- valued for each values of ^13. At the same time, 
we found that the smallest mass eigenvalue is rather stable for the allowed ^^13 in both cases 
(normal or inverted mass spectrum). We thus get absolute lower bound on the smallest 
mass as follows. 

• the case of normal mass spectrum: 

mi > 1 X 10"^ eV (3.15) 



12 




FIG. 1: A schematic view of neutrino mass spectra to show the obtained results. The horizon- 
tal coordinate corresponds to mi and the vertical coordinate corresponds to m2 and m^. The 
solid curved lines show the constraints to mass-squared difference obtained by neutrino oscillation 
experiments and the crossing points with the vertical dotted lines show the obtained results in 
Eqs.(jS31),(|HISl),(Eini) and (l8Tni . (l3T2|l . (ITOl) . 

• the case of inverted mass spectrum: 

ms > 1 X 10"^ eV (3.16) 

The lower bound for the inverted mass spectrum is very encouraging for the search for 
neutrinoless double beta decay j^]. We also have found that for certain value of 6*13 the 
Wilkinson microwave anisotropy probe (WMAP) results J2i^i < 0.70 eV, is violated 
in both cases. 



IV. RG EVOLUTION OF DIRAC NEUTRINO MASS MATRIX AND REANAL- 
YSIS OF THE NEUTRINO MASSES 

In previous sections, we derived the absolute values of neutrino masses, assuming JUjj 
coincide with up-type quark mass matrix at the weak scale. Such coincidence, however, 
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is realized only at the GUT scale, and the relation between Ain and up-type quark mass 
matrix is modified by the renormalization group (RG) effect. In this section, we thus improve 
the analysis of previous sections by including the RG effect. In the computation of RG 
evolution of Aio, we pay special attention to the successive decoupling of heavy right- 
handed neutrinos. 

In the following analyses, the framework we work on is a supersymmetric SO (10) GUT. 
We use a chain of spontaneous gauge symmetry breaking in which 5*0 (10) GUT breaks into 
the standard model via the minimal supersymmetry standard model (MSSM) : 

SO{10) ^'-^^ SU{5) ^'-^^ SU{3) X SU{2) x f/(l) ^'-^^ SU{3) x f/(l), (4.1) 

and for simplicity we assume that the symmetry breaking scales of 50(10) and SU{5) are 
sufficiently close to each other, i.e. (16//) ^ (45//) Jsi^ . We need to introduce multiple 
number of 10-dimensional repr. of Higgs field, since otherwise the up-sector and down- 
sector of quarks and leptons have degenerate masses (at GUT scale). Let us note that such 
problem is automatically evaded in our scheme, as supersymmetry inevitably demands two 
10-repr. of Higgs fields. The supersymmetry, of course, plays an important role in the 
unification of gauge couplings at the GUT scale. 

At first, we evaluate the evolutions of Yukawa coupling matrices of up- and down-type 
quark sectors from the weak scale towards the GUT scale (up-stream) in MSSM scheme 
and identify the obtained Yukawa couplings at Mqut with those of neutrino and charged- 
lepton sectors at the GUT scale, respectively. Then, we evaluate the evolutions of the 
Yukawa coupling matrices of lepton sectors from the GUT scale towards the weak scale 
(down-stream), taking into account the successive decoupling of three singlet neutrinos at 
their threshold energies. Using the obtained M.^ (the neutrino Yukawa coupling matrix 
X V ■ sin/5/v^), we can immediately write down the effective neutrino mass matrix at low 
energy. This mass matrix, however, still have three parameters, i.e. 0125 0235 013, after we 
fix R and to carry out the RG evolutions. To fix these three parameters, we use the 
conditions for three mixing angles, Eqs. ()2. lOj) . ()2. 1 1|) . ()2. 12|) . We then just follows the method 
we developed in the previous sections to get the eigenvalues oi M.y and compare the obtained 
ratio of mass-squared differences with the observed one shown in Eq. ()2.1Hj) . 

In MSSM, the RG equations of the Yukawa coupling matrices at one-loop are given by 
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dt 



IGvr 
In (/i/M,) , 



(a = u, d, u, e) 



(4.2) 



T„ = = 3 Tr[YuY^] + Tr[YX] = = 3 Tr[YdYj] + Tr[YeY^^] 
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Note that the coefficient we use is different from that used for Dirac neutrinos. This is 
because in the seesaw mechanism the presence of the inverse right-handed Majorana mass 
matrix modifies the coefficient as = 3 — 1. One problem here is that the evolutions 
of the Yukawa coupling matrices of quark and lepton sectors are coupled with each other 
since the term T^, the contribution of Higgs self-energy diagram, contains both of quark and 
lepton Yukawa couplings. It seems that when we solve the "up-stream" RG equation for 
quark Yukawa couplings, we need to know the values of leptonic Yukawa couplings, which 
are unknown. Fortunately, it turns out that in the estimation oi M.^ at the weak scale 
we can safely neglect the interference from another sector, as we show below. Namely, the 
Eq. ()4.2j) formally can be solved as 



Ya{t) = exp 



YaiO). 



(4.3) 



Since the exponent is of 0{a), we can expand this solution perturbatively up to its leading- 
order as 



Ya{0). 



Applying this arguments to the Yukawa coupling matrices of up-type quark and neutrino 
sector, and identifying these matrices at the GUT scale (/x ~ 2 x 10^^ GeV) , i.e. Yui^tGur) = 
YuitGur), Yi,{0), of our interest, can be expressed as 



n(o) 



((r.-G.)i + E^'^6n^) n(w) 

_i-^f^^^t((T.-G.)i + E^^nrj)_ 
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1 + /J"" - 1 + Y.ict - ci) nn^' 



b 

where we have used the relation Ty = T^- Therefore, the problematic factor has disap- 
peared, and we can treat the RG equations for quark sector and lepton sector independently. 
Adopting the input data set of the Yukawa coupling matrices of quark sectors at the weak 



scale (/i ~ Mz) 



1- v-sinf] ,. , 

mi = diag[Yu)ii x — -^=— [i = u, c, t) 



V ■ cos (3 

rriu = 2.33 MeV m, = 677 MeV rrit = 181 GeV 



ruj = diag(Yd)jj x ^= — (j = d, s, h) 



TUd = 4.69 MeV m, = 93.4 MeV = 3 GeV 

sin 012 = 0.2205 sin 623 = 0.0373 sin ^13 = 0.003 

, we can carry out solving the up- and down-stream RG equations. 

In order to know the threshold energies of decoupled singlet neutrinos, we need to know 
the values of R and m^. We scan the following ranges of these parameters, 

1 X 10"^ < i? < 1 X 10^^ (4.4) 
1 X 10^^ < mf < 2 X 10^^ (GeV), (4.5) 

and set tan/5 = 10 as a reference value (in fact, the obtained values are stable for tan/5 = 
1, 10 and 30). There are no difficulty in solving the up-stream RG equation. In the case of 
down-stream RG equation, however, we have to take into account the decoupling effects of 
singlet neutrinos at their threshold energies mf , ■ R and ■ R^. 

To clarify our prescription, we write the neutrino Yukawa coupling matrix as 

3 

UjR(j) . (4.6) 

In the energy region, mf < fi, the evolution of Yu is calculated by use of full members of Y^. 
Once /i passes through m^, i.e. in the region of ■ R < ^ < m^, the evolutions of Y^, is 
calculated by use of the subset of Y^ obtained by removing 3rd column of Y^, in the product 
of Y^,, i.e. {YyYj)ij = {Yy)ii{Y*)ji + {Yy)i2{Y*)j2, since u^r is decoupled in this region. We 
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can repeat the same procedure at each step of decouphng of single neutrinos, and we use 
the obtained at Mw to calculate J^u- 

2 2 

The obtained results are shown in FigI21and FigI21 where we define A = "^l . All of 

them satisfy the condition that m2 < ma, i.e. normal mass spectrum. According to these 

figures, A is rather insensitive to the values of m^. On the other hand, it is interesting to 

know that the allowed region for R is restricted to the range around 1 x 10~^. Since we 

fix the effective neutrino mass matrix at low energy completely, we can also analyze the 

absolute values of neutrino masses simultaneously and the obtained results are shown in 

FiglU We can find in this figure that the absolute values of lightest neutrino mass are stable 

significantly around certain region, 5 x 10^'^ < mi < 1 x 10~^ (eV). Note that these are 

higher than what we obtained in previous section (Eg. ()H.4|l ). 

m?(GeV) 

—■—I — 1 1 — 1— I — 1 1 — 1— I — 1 1 — 1— I — 

X X - 

X X - 

X 

X X 

X 

X 

X X 
« X - 

X X 
X X 
X X - 

10'^ 10"^ lO""^ 10'^ 

R 

FIG. 2: The horizontal coordinate corresponds to K and the vertical one corresponds to m^. 

2 2 

At the lattice points marked by cross, the obtained values of A = — | ^ is compatible with the 

allowed range 7.9/280 < A < 8.5/180 {3a C.L.) coming from neutrino oscillation experiments. 



2x10 
1x10 



1x10 
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1x10 



14 



1x10 



13 



17 



1 



10" 



i' ! ii'i 



I,, I II I" ' I ,, ! ' 



10" 



10"^ 10"^ 10""^ 10"^ 

R 

FIG. 3: The obtained values of A = as the function of R. Multiple points for each 

R correspond to different values of m^. We have taken tan/3 = 10. The horizontal solid line 
corresponds to the best fit value of A from experiments, A = 8.2/220, and the horizontal dotted 
lines show the allowed range of A at 3a C.L., i.e. 7.9/280 < A < 8.5/180. 

V. CONCLUSION 

In this paper, we derive the absolute values of three neutrino masses only invoking to the 
seesaw mechanism collaborated by the unification of neutrino Dirac mass matrix with that of 
up-type quarks. Especially for ^13 ~ 0, the obtained results are mi : m2 : rris a; 1 : 3 : 17 or 
1 : 50 : 250 in the case of normal mass spectrum of neutrino masses and mi ~ m2 : ms 2 : 1 
or 12 : 1 in the case of inverted mass spectrum and we found that there exists the lower 
bound of the smallest neutrino mass for the allowed range of ^13 at 3a C.L. In addition, 
we reanalyze the neutrino Yukawa coupling matrix at low energy including the effect of 
the RG evolution, paying a special attention to the successive decoupling effects of three 
singlet neutrinos and derive the ratio of mass-squared differences applying a method derived 
in previous sections. In this analyses, we find that in certain region of R there exist the 
solutions which are compatible with the experimental data of neutrino oscillations and that 
the absolute values of neutrino masses are higher that what we obtained in previous sections. 
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mi (eV) 



0.014 
0.012 
0.01 
0.008 
0.006 
0.004 
0.002 




10" 



10 



-5 



10 



-4 



10" 



FIG. 4: The obtained values of mi as the function of R. The horizontal coordinate corresponds 
to R and the vertical coordinate corresponds to mi (eV). The points marked by cross satisfy the 
condition, 7.9/280 < A < 8.5/180 {3a C.L.), and the other do not satisfy this condition. 

In the process, we impose some suitable assumptions for simplicity because our main 
purpose is to clarify the procedure to derive the neutrino masses, where large mixing angles 
observed in the leptonic sector is achieved invoking to the unification inspired by GUTs. 
Let us recall that, for instance, in 5*0(10) GUT, which incorporate left-right symmetry 
SU{A)ps X SU{2)l X SU{2){i, the neutrino Dirac mass matrix is naturally unified with that 
of up-type quarks and they are identical at least at the scale of Mqut- 

However, there still remain some questions about the justification of these assumptions. 
For example, it is an interesting question to ask, though difficult to answer, how the obtained 
results are modified when we take into account the 6 CP violation phases which in general 
embodied in Majorana mass matrix as physical degrees of freedom resulting from seesaw 
mechanism. 

These remaining questions will be addressed in ref.js^. 
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In general, it is also not so easy to solve non-linear coupled equations numerically. 
In the Yukawa sector, vev of heavy Higgs, (126//) = ((10,1,3)) ((10,1,3) isarepr. of5C/(4)p5X 
SU (2) L X SU (2)r), gives a Majorana mass to ur. In addition, vevs of two light Higgs, (10//) = 
((1, 2, 2)) (this is a 5'[/(4) singlet), give the Dirac masses to quark and lepton sectors and we 
can achieve the desirable relations, niu = rriy and rrid = nie at the GUT scale 
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